Anomalous scaling dimensions and stable charged fixed-point of type-II 

superconductors 
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The critical properties of a type-II superconductor model 
are investigated using a dual vortex representation. Comput- 
ing the propagators of gauge field A and dual gauge field h in 
terms of a vortex correlation function, we obtain the values 
77a = f and 77h = f for their anomalous dimensions. This pro- 
vides support for a dual description of the Ginzburg-Landau 
theory of type-II superconductors in the continuum limit, as 
well as for the existence of a stable charged fixed point of the 
theory, not in the ZDXY universality class. 

PACS numbers: 74.6Q.-w, 74.2Q.De, 74.25.Dw 



Determining the universality class of the phase- 
transition in a system of a charged scalar field coupled 
to a massless gauge field, such as a tyne-II superconduc- 
tor, has been a long-standing problemEI. Analytical and 
numerical efforts have recently focused on the use of a 
dual description of the Ginzburg-Landau theory|,(GLT) 
of type-II superconductors, pioneered by KleinertB, in in- 
vestigating the character of a proposed novel stable fixed 
point of the theory for a charged superconducting con- 
densate, in which case the 'iDXY fijed point of the neu- 
tral superfluid is rendered unstabletl The dual formu- 
lation has also been employed to investigate the possibil- 
ity of novel broken symmetries in the vortex liquid phase 
of such systems in magnetic field^'El. 

The GLT is defined by a complex matter field ij) cou- 
pled to a massless fluctuating gauge field A with a Hamil- 
tonian 



H^^A. — 



2(Vx 



(1) 



Here, e is the electron charge, and i?^,A is invariant 
under the local gauge-transformation -0 — > ■0exp(i0), 
A ^ A + V0/2ie. The GLT sustains stable topological 
objects in the form of vortex lines and vortex loops, the 
latter are the critical fluctuations of the theoryoB. These 
objects are highly nonlocal in terms of but a dual 
formulation offers a local field theory for them. The con- 
tinuum dual representation of the topological excitations, 
{in D — i only), consists of a camplex matter field cou- 
pled to a massive gauge field hn, with coupling constant 
given by the dual charge ed, and with dual Hamiltonian 

T 



i(Vxh)^ 



\ (V X A)^ 



ie (V X h) • A. (2) 



The massiveness of h reduces the symmetry to a global 
C/(l)-invariance. For details on how to obtain this dual 
Hamiltonian, we refer the reader to the thoroughj^xpo- 
sition of this presented in the textbook of KleinertQ. For 
e ^ the original GLT in Eq. Jl| has a local gauge sym- 
metry, the dual theory in Eq. g has a global U{1) sym- 
metry. In the limit e ^ 0, A decouples from ?/; in Eq. [|, 
describes a neutral superfluid, and the symmetry is 
reduced to global C/(l). The dual Hamiltonian i/^.h de- 
scribes a charged superfluid coupled to a massless gauge 
field h with coupling constant e^, and the global sym- 
metry is extended to a local gauge symmetry. Hence, 
when e ^ 0, the dual of a neutral superfluid is isomor- 
phic to a superconductor. Integrating out the A field 
in Eq. ^ produces a mass-term e^h^/2, where an exact 
renormalization-atoup equation for the mass of h is given 
by de^/dlnl = e^B. Therefore, when e ^ 0, then — > oo 
as ^ ^ oo. This supresses the dual gauge field, and the 
resulting dual theory is a pure |(/)|'*-theory. Hence, in 
the long-wavelength limit, the dualjof a superconductor 
is isomorphic to a neutral superfluiaa. 

In this paper, we obtain the |-anomalous scaling di- 
mensions r]A of the gauge fieldO'EI, as well as r]^ of the 
dual gauge field, not previously considered, directly from 
large-scale Monte-Carlo simulations. At a 3DXY critical 
point, 77a = ?7h = 0. We find that {rjA = l,77h = 0) when 
e ^ 0, and that (77A ~ 0,ilh — 1), when e ~ 0. We also 
contrast the anomalous dimension of the dual mass field 
(j) at the dual charged (original neutral) and dual neutral 
(original charged) fixed points, obtaining 77^ ~ —0.24 in 
the former case, and 770 = 0.04 in the latter. 

A duality transformation, to a set of interacting vortex 
loops, is performed on the London/ Villain approximation 
to the GLT. In this approximation the partition function 
is 



■' {n} L X 1^ 



^ (A^ - eA - 2™)^ 



(3) 



l(V-iedh) 



Here, 9 is the local phase of the superconducting order 
parameter while n is an integer- valued velocity field 
(not vortex field) introduced to make the Villain poten- 
tial 27r-periodic. The symbol A denotes a lattice deriva- 
tive. Amplitude fluctuations are neglected in this ap- 
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proach. The validity of this approximation for 3-D sys- 
tems, has recently been iiattestigated in detail, both nu- 
merically and analyticallyty. 

An auxiliary velocity field v linearises the kinetic en- 
ergy. Performing the ^-integration constrains v to sat- 
isfy the condition A • v = 0, explicitly solved by writing 
V = A X h, where h is forced to integer values by the 
summation over n. Introducing an integer-valued vortex 
field m = A X n, and using Poisson's summation formula, 
we find 

S{A, h, m) = ^ Unim • h + ^ (A x h)^ 

-Hie(Axh)A+i(Ax A)^|. (4) 

Integrating the gauge field in Eq. ^ produces a mass term 
giving an effective theory containing the vortex 
field m coupled to a massive gauge field h 

Z{(3,e) = J Dh^Y[SA-mfiexp -^|27^^m•h 



(AqA_q) = 



Q|- + m2 



47r2/3mgG'(q) 
IQl'flQl'+m2 



2/3 



Ql' + m2 



\Q\' + ml 



(8) 



(9) 



where G(q) (mqm_ 



A ^ and 



I - 

g-jq-fi correlation functions have been calculated in 
the transverse gauge V • A = V • h = 0. Both of the fields 
h and A are renormalized by vortex fluctuations, albeit 
in quite different ways. 

Invoking the standard form (q^ + m^^) for the cor- 
relation functions in the immediate vicinity of the critical 
point in the limit g — s- 0, we find the following expressions 
for the effective masses, 



m, 



effj 



lim 



9-0 1 + 47r2/3G(q)g-2 ' 



{m} X 



1 + 



9-0 2^ (^^ _ 2.^^G(q)^ 



(10) 
(11) 



(5) 



When e 7^ the correlation function for A assumes the 
form 



The variables m in Eq. || describe a set of interacting 
vortices, where the interactions are mediated through the 
gauge field h. The variables in Eq. || are defined on 
a lattice which is dual to the lattice from Eq. |3[ and 
the behavior with respect to temperature is inverted in 
the new variables. The 6 field in Eq. ^ describes order, 
while the m field represents the topological excitations 
of the 9 field. These excitations destro¥,superconducting 
coherence, and hence quantify disordem. 

Integrating out the h field in Eq. ||, we obtain the 
Hamiltonian employed in the present simulations, 



(AqA_q) 



1 



q- 



(12) 



at the critical point. To determine 77A, we compute the 
vortex correlator G{q). For A << L = 40 , we expect the 
following behaviour for G(q) in the limit g — > 0, 



T<T,^ G(q) cx q^, 
r-T,^G(q)(x g", 
r>r,^G(q)oc C{T). 



(13) 
(14) 
(15) 



When these limiting forms are inserted in Eq. nO, we see 



H{m) = -271 Jo m(xi)V^(xi - X2)m(x2), (6) ^^^^t for T < T^, will be finite through the Higgs 



Xl,X2 

g-iqx 



r4E^sin2(f)+A- 



(7) 



In Eq. ^ the charge e and lattice-spacing a have both 
been set to unity, and A is the bare London penetration 
depth. At every MC step, we attempt to insert a loop 
of unit vorticity and random orientation. A new energy 
is calculated from Eq. 6, and the proposed move is ac- 
cepted or rejected according to the Metropolis algorithm. 
This procedure ensures that the vortex lines of the syj 
tem always form closed loops of random size and shap' 
In all simulations, a system size of 40 x 40 x 40 was used, 
and up to 1.5 TO^ sweeps over the lattice per temperature 
were used. 

To investigate the properties of A and h at the charged 
critical point of the original theory, Eq. pi we have cal- 



cff 

Mechanism (Meissner effect). For T > Tc we will have 
= as in the normal case of a massless photon. 
Assuming G{q) cx q^ precisely at the critical point, it is 
seen that 7] corresponds to 7/a from Eq. |l2|. We thus 
identify the scaling power of G(q) at the critical point 
with the anomalous dimension of the massless gauge field 



culated the correlation functions (AqA_q) 
in terms of vortex correlations, obtaining 



and 



;hqh_ 
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FIG. 1. log-log plot of G{q) for the three alternatives in 
Eq. |l| - with A = a/2. For this A, Tc = 1.446. Apart 
from the point q = qmin,r ~ 1.446 the error bars are smaller 
than the symbols used. 



All three limiting forms Eqs. [13[ - |lq are shown in Fig. 
l|. The gauge field masses m^g and to^ in Eqs. ^ and 
are shown in Fig. ||. At the critical point G{q) cx q, so 
that rjA. — 1. Note that, while vanishes at T = Tc, 
rn^ff is finite but non-analytic. As a result of the vortex 
loop blowout, the screening properties of the vortices are 
dramatically increased, and rri^g increases sharply. 




FIG. 2. m^/mo and m^ff/mo as functions of T. 

To find 77h independently, we consider first the un- 
charged case A —^ oo, mo 0. First, at an intermediate 
step in the transformation Eqs. ^ - ^, the action reads 

S{(3, e) = - X] I + e»A . 1 + i (V X A)^! . (16) 
Here, 1 is an integer field of closed current loops. Setting 



e = in Eq. the action of the dual Villain model is 
obtained, 

~Sv{P, r) = |27rzm • h + ^ (A x h) V ^m^j. (17) 

Here, a term Fm^/2 has been added, and S'y(/3, F) cor- 
responds to the Villain- action in the limit F ^ 0. How- 
ever, it is physically reasonable to propose that the limit 
F — > is non-singular, since the added term is short- 
ranged. It should therefore be an irrelevant perturba- 
tion, in renormalization group sense, to the long-ranged 
Biot-Savart interaction governing the fixed point, which 
is mediated by h. Rescaling h he/27r in Eq. 17, we 
have0Z(/3,e) = Zy (eV47r2, 1/2^, leaving Eqs. ||and 
p7| interchangeable; ?7h from Eq. |l^ should have the same 
value as r/A from Eq. |l^. The above is demonstrated by 
our simulations based on Eqs. |^^, which are independent 
of the proposed form Eq. |7^. 

To determine Jjh we study the correlation function 
(hqh_q) (Eq. ^ in the limit toq 0. At the un- 
charged fixed point of the original theory, which is 
the charged fixed point of the dual theory, we have 
lim,^0 27r/32G(q) = {l-C2{T))q' + ...,q' ^C:,{T)q'+^^ + 
and g2 _ Ci{T)q* + for T < T^,T ^ T^, and 
T > Tc, respectively. Here, C^iT) corresponds to the 
helicity modulus (superfiuid density )t3, C^iT) is a crit- 
ical amplitude, and Ci{T) is the inverse of the mass of 
the dual gauge field for T > Tc. Correspondingly, we 
have lim,^o < hqh_q >= 2pC2/q^, 2(30:!./ q^'"^^ , and 
2[3Ci, for T < Tc,T = Tc, and T > Tc, respectively. 
Note that h is massless for T < Tc, while it is massive for 
T > Tc, the dual system exhibits a "dual Meissner-effect" 
for T>Tc. At T = Tc, we have q^ (hqh_a) ~ C3{T)q''^. 
A plot of q^ (hqh_q) is shown in Fig. y. A linear be- 
haviour at T — Tc is found, implying that jyn = 1 when 
e = 0. Since r^h = 1 in the uncharged case, this provides 
further support for the Hamiltonian Eq. |^. 

We now set e ^ 0. The gauge field h becomes massive 
via the term e^h^/2, which appears after integrating out 



the A field in Eq. 



In this case, lim„^o(hqh_ 



2/3/m§ from Eq. ^, and h(r) would naively have the 
trivial scaling dimension (2 — d) /2. However, the mass 
term offers us a freedom in assigning dimensions to e 
and h, by introducing renormalization Z-factors, here 
e' = Zl/^e and h' = Z~^^\. 

Prior to integrating out A in Eq. |^, the mass ap- 
pears in the term ie (V x h) • A. Integration of the (f> 
field, partial or complete, can only produce (V/i — Cdh)- 
terms. In particular, this must hold during integration 
of fast Fourier- modes of the cj) field. Thus, the term 
i(V x h) • A is renormalisation group invariant, i.e. its 
prefactor must be dimensionless. In terms of scaled fields, 
at the charged fixed point of the original theory, we have 
A' = ^A^^^A, with Za oc Z''^, 77a = li. For h, we use 
Zh I , where A is not an anomalous scaling dimen- 
sion (h is massive, cf. Fig. H), but rather a contribution 



3 



to the engineering dimension of h. Inserting this into 
the crossterm ie (V x h) • A, we find the scaUng dimen- 
sion (77a + ^) /2 — 1, which must vanish. This gives the 
constraint A = 1 to avoid conflicting results for rjA.. 

Remarkably, therefore, the scaling dimension of h at 
T = Tc is the same in both cases tuq — and mg =/= 0. 
The results for 77A and rjh in the previous paragraphs, 
are summed up in Table 0. 
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FIG. 3. g^(hqh_q) for two different T. For A = 00, 
Tc = 3.00. 

We next consider the distribution of vortex loop sizes in 
the model Eq. connecting the vortex loop distribution 
to the anomalous dimension of (j) a-t Tc both for the case 
e = and e ^ Q. During the simulations, we sample the 
distribution of loop-sizes D{p), where p is the perimeter 
of a lUwp- This distribution function can be fitted to the 
formO'D 

D{p) (xp-"e~^P''^'^\ (18) 

where s{T) is an effective line-tension for the loops. Fig- 
ures showing the qualitative features of D{p) can be 
found in Ref. |[ The critical point is characterised by 
a vanishing line-tension, and close to the critical point 
we find that e{T) vanishes as £(T) cx |T - Td^*. 

The vortex loops are the topological excitations of the 
GL and 3DXY models, at the same time they arc the 
real-space representation of the Feynman diagrams of the 
dual field theory. By sampling D{p), we obtain infor- 
mation about the dual field 0, particularly 7^ can he 
identified as a susceptibility exponent for the cj) fieldQ. 
Using the scaling relation 7^ — i/^{2 — ry^), and the im- 
portant observatioijuthat even at the charged dual fixed 
point Vfj, = PsDXYa, this also gives us a value for the 
anomalous scaliag dimension rj^ when we use the value 

1^3DXY = O.67M. 

In Ref. H the vortex loops of the 3DXY model have 
been studied meticulously, yielding the value ti^(0) = 
—0.18 ± 0.07. Since the dual of this model is isomorphic 



to a superconductor, r],p{0) should be similar to r]^{e) of 
the original GLT. 
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FIG. 4. Ine(r) as a function of In \T - Tc]. The upper line 
shows the charged case with finite e, and the lower line shows 
the neutral case with e — 0. The slopes of the two straight 
lines are 7^ = 1.315 and 7^ = 1.51, corresponding to the 
anomalous dimensions rj^ = —0.24 (neutral, i.e. dual charged) 
and ri^ = 0.04 (charged, i.e. dual neutral), respectively. 

We have studied the vortex loop distribution in both 
the neutral and the charged case. In the former case 
we find rj^ ~ —0.24, in good agreement with Ref. ||. In 
the latter case the dual theory has a U{1) symmetry, 
and we would expect to find 77^ — rj^jjxY- The exponent 
113DXY has recentljf-been determined with great accuracy 
to r]3DXY = 0.038t^, whereas we findn^ ~ 0.04 which 
compares well with this value. Fig. H shows e(T) for 
both the charged and uncharged models. It is evident 
that they belong to two different universality classes. 

In the case e ^ 0, which corresponds to the dual neu- 
tral case, the inverse 0-propagator is given by = q'^ + 
E(g), where E is a self-energy, and S(g) ~ g^"'' by defi- 
nition. This gives a leading order behavior G ^ 
provided rj > 0, and we find 77 = 0.04 for this case. On 
the other hand, for the case e = 0, which corresponds 
to the dual charged case, dual gauge field fluctuations 
alter the physics, softening the long- wavelength (j> field 
fluctuations. We obtain G^^ = + ^(g), again with 
I](g) ~ which now gives a leading order behavior 

G ~ 1/q-^^^, provided 77 > —2. Our result 77 = —0.24 for 
the case e = (dual charged) is consistent with this, and 
also with the absolute bounds 77>2 — _D = — l,inr' = 3. 

A consequence of the above is that in Z? = 3 dimen- 
sions, A ^(c-2)/(2-r,A) ^ ^ at the charged critical 
point, in contrast to A ~ at the 3DXY neutral crit- 
ical point. Since our results have been obtained directly 
by MC simulations, they are valid beyond all orders in 
perturbation theory. 
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FP, original theory 




FP, dual theory 








Neutral 3DXY 


1 


Charged 


Finite 


1 


Charged 





Neutral 3DXY 



TABLE I. Values of r]A and r/h at the stable neutral and 
charged critical points of the original and dual theories. FP 
is an abbreviation for fixed point. 
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